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Abstract. We discuss the time-convolutionless (TCL) projection operator approach to transport in closed 
quantum systems. The projection onto local densities of quantities such as energy, magnetization, particle 
number, etc. yields the reduced dynamics of the respective quantities in terms of a systematic perturbation 
expansion. In particular, the lowest order contribution of this expansion is used as a strategy for the 
analysis of transport in "modular" quantum systems corresponding to quasi one-dimensional structures 
which consist of identical or similar many-level subunits. Such modular quantum systems are demonstrated 
to represent many physical situations and several examples of complex single-particle models are analyzed 
in detail. For these quantum systems lowest order TCL is shown to represent an efficient tool which also 
allows to investigate the dependence of transport on the considered length scale. To estimate the range 
of validity of the obtained equations of motion we extend the standard projection to include additional 
degrees of freedom which model non-Markovian effects of higher orders. 



PACS. 05. 60. Gg Quantum transport 
and irreversible thermodynamics 



05.30.-d Quantum statistical mechanics - 05.70.Ln Nonequilibrium 



1 Introduction 

■ In recent years the application of projection operator tech- 
niques [1I2I3I4I5I6] to transport investigations in (closed) 
quantum systems has been suggested |7I8I9] . Within these 
approaches a suitable projection onto local densities of 
pertinent transport quantities is used in order to obtain 

■ the reduced dynamics of these quantities in terms of a 
systematic perturbation expansion, typically w.r.t. some 
(small) interaction strength. For certain examples the low- 
est order contribution of this expansion has been shown 
to yield reliable predictions on transport [8] as well as on 
its length scale dependence [7]. 

However, the lowest order truncation may become ques- 
tionable for strongly interacting quantum systems, of course 
But already in the weak coupling case the validity of this 
truncation generally is restricted to short time scales [TU] . 
This fact seems to be problematic, especially since the 
relevant time scale can be very long for small interaction 
strengths or large length scales [7], e. g., in the thermody- 
namic limit. For this reason the additional consideration 
of higher order contributions appears to be indispensable 
[5], contrary to the statements in, e. g., Ref. [9]. Because 
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such a consideration has not been provided in the litera- 
ture so far, the main intention of the present paper is the 
incorporation of higher order terms as well, particularly 
the next-to-lowest order contribution. 

To start with we will introduce the notion of a modu- 
lar quantum system in Sec. [2] Modular quantum systems 
are demonstrated to represent many physical situations 
and several examples will be given. In the context of these 
quantum systems the time-convolutionless (TCL) projec- 
tion operator technique [31514) is subsequently discussed. 
In Sec. [3] the projection onto local densities and lowest or- 
der TCL is firstly shown as an appropriate method which 
also allows to investigate the dependence of transport on 
the considered length scale. In particular explicit condi- 
tions for the applicability of the introduced method are 
given. The next Sec. [5] is concerned with the higher order 
contributions of the TCL expansion and a suitable estima- 
tion is derived for the special case of interactions with van 
Hove structure |llll2U0j . In Sec. [5] this estimation is used 
in order to determine the range of validity of lowest order 
TCL and an interpretation in the context of length scales 
is provided. Section [5] finally applies the concepts of the 
previous Sections to complex single-particle models, e. g., 
to models with disorder. Our results are confirmed by nu- 
merical solution of the full time-dependent Schrodinger 
equation. 
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2 Modular Quantum Systems and Diffusive 
Dynamics 

In the present paper we consider so-called "modular" quan- 
tum systems which have a quasi one-dimensional structure 
and consist of N identical or at least similar many-level 
subunits. These subunits are described by a local Hamil- 
tonian /i u and the next-neighbor interaction between two 
adjacent subunits is denoted by Aw^+i, where A ad- 
justs the overall coupling strength. The total Hamiltonian 
H = H + A V is given by 



Ha 



N-l 

E v 

(i=0 



V 



N-l 

E 



fl,fi+l j 



(1) 



where we employ periodic boundary conditions, i. e., we 
identify /i = N with \i = 0. Such a description obviously 
applies to one-dimensional structures such as chains of 
atoms, molecules, quantum dots, etc. But also s = 1/2 
spin chains fit into this scheme of description, because 
a segment of the chain, that is, a number of spins and 
their mutual interactions can be chosen in order to form 
a suitable many-level subunit. Similarly, spin lattices or 
higher-dimensional models of the Hubbard type can be 
treated by the use of fl}, if a whole chain (2D) or layer 
(3D) is considered as a single subunit. Thus, for a large 
class of quantum systems an adequate way of description 
is offered by the Hamiltonian of Eq. (fl}, see Sec. [5] for 
details. 

In this paper we primarily deal with "single-particle" 
quantum systems, that is, those quantum systems which 
allow to restrict the investigation to a linearly instead of 
an exponentially increasing Hilbert space. This is mainly 
done in order to open the possibility for a comparison of 
the theoretical predictions with the numerical solution of 
the full time-dependent Schrodinger equation. Hence, we 
may suppose 



n-l 
i=0 



l/M) {n, 



(2) 



1 |^,i)( M +l,j|+H.c, (3) 



where n denotes the number of levels within a subunit. 
Without loss of generality, we have additionally assumed 
an off-diagonal block structure of the interaction, that is, 
i) and \n,j) are not coupled. 

Of particular interest is the local density p^ (t) of, e. g., 
energy, magnetization, excitations, particles, probability, 
etc. This quantity is expressed as the expectation value of 
a corresponding operator Pfl , 

Pn(t) = Tr{/?(t)j5 u } , = E^m I/MX/Ml > ( 4 ) 

where p(t) is the full system's density matrix. Since the 
operators p^ are assumed to be diagonal in the energy 



representation of the uncoupled system, we restrict our- 
selves to those quantities p^ (t) which are conserved for the 
special case of A = 0. However, this restriction still allows 
to investigate transport for a large class of systems, as 
will be demonstrated in Sec. [5] Our aim is to analyze the 
dynamical behavior of the p^it) and to develop explicit 
criteria which enable a clear distinction between diffusive 
and other available types of transport, e. g., ballistic or 
insulating behavior. 

The dynamical behavior may be called diffusive if the 
p^t) fulfill a discrete diffusion equation 



p fl (t) = D[p^ 1 (t)-2p fl (t)+p fl+1 (t)\ 



(5) 



with some \x- and t- independent diffusion constant D. It is 
straightforward to show [multiplying ([5]) by respectively 
/^ 2 , summarizing over /i and manipulating indices on the 
r.h.s.] that the spatial variance 

N-l [N-l 1 2 

a 2 (t) = Y / V 2 Pn(t)- £>P/.(*) (6) 

fj,=0 L "=0 

increases linearly with t, namely a 2 (t) = 2Dt. By con- 
trast, ballistic behavior is characterized by cr 2 {t) oc t 2 , 
whereas insulating behavior corresponds to <J 2 (t) = const. 

According to Fourier's work, diffusions equations are 
routinely decoupled with respect to, e. g., cosine-shape 
spatial density profiles 



N-l 



Pq(t) = C q E C0S ('?A i )PA'( < ) J 9 
u=0 



2?rfc 



(7) 



with k — 0, 1, ... , A^/2 and a suitable normalization con- 
stant C q . Consequently, Eq. © yields 



p g (t) = -2(1- cos q)Dp g (t). 



(8) 



Therefore, if the quantum model indeed shows diffusive 
transport, all modes p q (t) have to relax exponentially. If, 
however, the p q (t) are found to relay exponentially only for 
some regime of q, the model is said to behave diffusively 
on the corresponding length scale I = 2ir/q. One might 
think of a length scale which is both large compared to 
some mean free path [below that ballistic behavior occurs, 
a 2 (t) tx t 2 ] and small compared to, say, some localization 
length [beyond that insulating behavior appears, o 2 (t) = 
const.], see Sec. [SI 



3 Projection onto Local Densities and Second 
Order TCL 

A strategy for the analysis of the dynamical behavior of 
the p q (t) is provided by the time-convolutionless (TCL) 
projection operator technique 3 5 4J. This technique, and 
the well-known Nakajima-Zwanzig (NZ) method as well, 
are applied in order to describe the reduced dynamics of 
a quantum system with a Hamiltonian of the form H = 
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Hq + XV, see Refs. [12. Even though projection operator 
methods are well-established approaches in the context 
of open systems, the following application to transport 
in closed systems is a novel concept. Let us remark that 
the application of these techniques only requires that the 
pertinent observables commute with Hq. 

Generally, the full dynamics of a quantum system is 
given by the Liouville-von Neumann equation 



d_ 



p(t)=-i[\V(t),p(t)]=£(t)p(t), 



(9) 



where time arguments refer to the interaction picture. In 
order to describe the reduced dynamics of the system, one 
has to construct a suitable projection operator V which 
projects onto the relevant part of the density matrix p{t). 
In particular, V has to satisfy the property V 2 = V . Be- 
cause in our case the relevant variables are the local den- 
sities p^{t), we choose 



N-l 



N-l 



This choice indeed fulfills the property V 2 
additionally normalize 

n-l 

Tr{(^) 2 } = EH) 2 = 1 > 



(10) 



■p, if we 



(11) 



which can be done without loss of generality. Note that 
this normalization typically implies p 1 ^ ~ 1/ \/n, at least 
for the models in Sec. El 

Once some projection operator has been defined, the 
TCL formalism routinely yields [not only for our choice 
of V in Eq. (fTO]) ] a closed and time-local equation for the 
dynamics of V p(t), 



d_ 

dt 



V p{t) = JC(t) V p(t) + J(t) (1 - V) p(0), 



(12) 



and avoids the often troublesome time-convolution which 
appears, e. g., in the context of the NZ method. Eq. (Tp2"j) 
and, in particular, its time-locality are a standard result 
and a direct consequence of the TCL formalism [31514] . 
The time-locality of the TCL equation may be understood 
as the result of a forward and backward propagation of 
the corresponding NZ equation in time. In Sec. [5] we will 
briefly comment on the NZ approach and, especially, on 
its implications for our class of models. 

For initial conditions p(0) with V p(0) = p(Q) the inho- 
mogeneity I{t) on the r.h.s. of (|12[) vanishes. [But for the 
models in Sec. [6] there are numerical indications that I{t) 
is even negligible for other p(0), see also Refs. |13I5I10| .] 
The generator K,(t) is given as a systematic perturbation 
expansion in powers of the coupling strength A, 



lC{t) = £ A m /C m (£). 



(13) 



The odd contributions of this expansion vanish for many 
models and for our models as well, that is, K,2m-i = 0. 



Consequently, the lowest non- vanishing contribution is the 
second order K, 2 {t) which reads 



ic 2 (t) 



f dt x V C(t) £(ti)P . 
Jo 



(14) 



The next non-vanishing contribution is the fourth order 
K-4(t) which is given by 



£ 4 (t): 



t rt± i<t2 

dh / di 2 / di 3 



-rc(t)c(t 2 )vc(h)c(t 3 )r 

-VC(t)C(t 3 )VC{h)C(t 2 )V 



(15) 



The truncation of the generator (fT3")) to lowest order, 
i. e., its approximation by the second order (|14p . is com- 
monly done in the case of small A but is in general re- 
stricted to short time scales. As already mentioned in 
the introduction, this truncation seems to be problem- 
atic, since the relevant time scale can be very long, es- 
pecially if A is small. However, the rest of this Section is 
only devoted to the second order truncation. In the fol- 
lowing Sec. |4] we will additionally discuss the fourth order 
(|15p and subsequently show that its incorporation is in- 
dispensable in order to estimate the accuracy of the sec- 
ond order prediction. Moreover, we will demonstrate that 
some non-diffusive transport phenomena such as localiza- 
tion can not be predicted correctly by a mere second order 
consideration, see Sec. H] 

Plugging the projector (JTOj) into ([12]) and ([14]) leads to 



(16) 



Note that the sum does not run over all v, because only 
adjacent subunits are coupled. The time-dependent rates 
R/i,v(t) are defined by 



(17) 



where we have introduced the correlation functions 



<Wr) = Tr{[p^V(t)][p v ,V(t 1 )]} , r = t-t x (18) 

with C^. v (t) — Cv.^t). [The trace is independent from 
the order of commutators and time arguments.] 

So far, Eq. ([Po]) is exact. For the further simplification 
of (|16p two assumptions have to be made now. 
(i.) We assume that 1 (i), Rn : n(t) and i? MiM+ i(t) de- 

pend only negligibly on the concrete choice of p. Hence, 
it follows that R n . n -i(t) w R u+lyU (t) = R u . u+1 (t). This 
is fulfilled exactly, if the system is translational invariant, 
e. g., if the coefficients in @, © and ^ are the same for 
each subunit. 
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(ii.) Moreover, we assume R^^-xit) w —R^^{t)/2. This 
is perfectly satisfied if the sum of all local densities is con- 
served since the overall conservation implies 

N-l 

Y J [P»,V{t 1 )]=0. (19) 

v=Q 

Consequently, the sum of C^ v {t) over all v, eventually 

2C w _i(r) + C^t), also vanishes. 

Due to (i.) and (ii.) the Fourier transform of (fT(3|) reads 

p q (t) = -WR(t)p q (t), W = 2(l-cosg)A 2 (20) 

with a single rate R(t) = —R fl!fl (t)/2. Note that this rate 
is still time-dependent. Remarkably, the dependence on q 
and A simply appears as an overall scaling factor W. 

The models in Sec. [5] typically feature a correlation 
function C(r) which decays completely within some time 
scale tc- After this correlation time C(r) approximately 
remains zero and R{t) takes on a constant value R, the 
area under the initial peak of C{r). Since the correlation 
time apparently is independent from q and A, it is always 
possible to realize a relaxation time tr oc 1/W which is 
much larger than tc, e. g., in an infinite system there 
definitely is a small enough q. For tc <C tr the second 
order truncation (|20p immediately yields 

p q (t) = -WRp q (t), (21) 

and the comparison with ([8]) clearly indicates diffusive be- 
havior with a diffusion constant D = A 2 R. 

However, the present method is not restricted to the 
investigation of diffusive transport phenomena and (|20[) 
is applicable as well in order to completely classify the 
dynamics of p q (t) which decay on relatively short time 
scales below tc or on very long time scales where, e. g., 
correlations may reappear, see Ref. [7] and Sec. El 

4 Extension of the Projection and Fourth 
Order TCL 

This Section and the next Sec. [5] as well are concerned 
with the question, why and to what extend the projection 
onto local densities and second order TCL yield reliable 
predictions for the dynamical behavior of the p q (t). The 
answer to this question surely requires the consideration of 
the higher order contributions of the TCL expansion, e. g., 
the investigation of the fourth order K,4(t), cf. (fT5|) . But 
already the direct evaluation of /C4 (t) turns out to be ex- 
tremely difficult in general, both analytically and numeri- 
cally. The problem is mainly caused by the first integrand 
in the expression (|15p which significantly contributes to 
A ;•:/ :. 

We therefore present an alternative approach which 
is essentially based on the following idea: The influence 
of the higher order terms may decrease substantially, if 
the projection additionally incorporates variables which 
are not of particular interest by themselves but poten- 
tially affect the dynamical behavior of the local densities. 



This idea is obviously only useful if the complexity of the 
higher order terms decreases to a larger extend than the 
complexity of the second order contribution increases. 

We concretely choose N (n — 1) additional variables 
c^Jt) which are also given as the expectation values of 
corresponding operators at. These operators are assumed 
to be diagonal in the energy representation of the uncou- 
pled system and are given by 

n-l 

al = J2^\^j)(^3\ (22) 

with fj, = 0, 1, . . . , N — 1 and i = 0, 1, . . . , n — 2. Addition- 
ally, the operators are supposed to fulfill 

Tr{a^} = <^-, Tr{&i l p li }=0 (23) 

such that the set of all a^, p M spans the whole space of 
diagonal matrices. The extended projector 

N-l n-2 

consequently projects onto the diagonal elements of the 
density matrix p(t). Thus, the complement 1 — V in the 
first integrand of (|15[) is a projection onto non-diagonal 
elements. But non-diagonal contributions are negligible if 
we restrict ourselves to interactions with the so-called van 
Hove structure, that is, V 2 is essentially a diagonal ma- 
trix, cf. Refs. [11|12|10] . (If the van Hove property is not 
fullfilled, the second order prediction is not reliable at all 
[10] •) The latter fact indicates that the extended projector 
transforms the largest part of original fourth order effects 
into the second order. For the general case of interactions 
without van Hove structure, however, the choice of diag- 
onal operators in Eq. (|22j) may not represent the most 
relevant variables for fourth order corrections, of course. 

Plugging the extended projector ([24]) into Jl2|), (fl"4f 
[and assuming (i.) translational invariance, (ii.) overall 
conservation of local densities] one finds 

p„(t) =A 2 R(t) [p M _i(t) - 2p M (t) +p„+i{t)} 

p+l n-2 

+ A 'E E (25) 

u—fi — 1 i—0 
fi+1 n-2 

W = a 2 £ + E ( 26 ) 

u=fj,-l j=0 

with time-dependent rates a l ^ v {t) and o^„(i) which are 
defined analogously to (fl7|) as integrals over correlation 
functions C % {t) and C^(r), of course. While C l {t) is 
given by 

Cl v (T)=Tr{[p^V(t)][&l, Hi)}}, (27) 

C % ^ v (t) is obtained by using a' and &(, as first arguments 
in the above commutators. For simplicity, however, we 



Robin Steinigeweg et al.: Projection Operator Approach to Transport in Complex Single- Particle Quantum Systems 5 



set the corresponding rates ofy^it) = for the following 
reason: Because we still consider initial conditions with 
a u(0) = 0, a significant increase of a l „{t) certainly arises 
only from the first part of (|26p. at least for sufficiently 
small times. 

The translational invariance implies that a AtjAl _i(t), 
and cv AtjAt +i(t) do not depend on the concrete choice 
of [i. (For clarity we suppress the fixed index i here.) 
But a MiM _i(£) w 
the additional assumption (iii.) a^^it) w u v ^{t), which 
is exactly fulfilled for mirror symmetry, that is, the co- 
efficients in © do not depend on the index order. The 
overall conservation of the local densities finally leads to 
cW-i(i) = _c w(*)/ 2 - 

Applying the Fourier transform to (|25[) yields 



Since the set of all a 1 , p forms a complete orthonormal 
basis, the fc-sum on the r.h.s. of the above equation is 
identical to Sij — p % pP . As a consequence the remaining 
sums over i and j can be performed independently from 
each other. Finally, by integrating over the independent 
variables t and n, a straightforward calculation leads to 



1 (t) only holds true if we make E «*(*) = E ~ R ® R ^ > ( 34 ) 



fe=0 



i=0 



p q (t) = -w\R(t) Pq {t) + «*(*) 4w 



i=0 



(28) 



and the Fourier transform which results only from the first 
part of (|26)) is given by 



(29) 



d q (t) = -Wa t (t)p q (t) 



with rates a l (t) = — alJt)/2. Finally, integrating (|29 
and inserting into ([28]) leads to 



(30) 



p g (i) = [-T^i?(i) + VF 2 S(<)] Pg (<), 



Remarkably, the extended projection has lead to an ad- 
ditional contribution with an overall scaling factor W 2 = 
4(1 — cos q) 2 A 4 . This simple scaling suggests the equiva- 
lence between small coupling strengths and large length 
scales. Especially, the factor A 4 indicates that the addi- 
tional contribution can be interpreted partially as a fourth 
order effect of the original projection. Note that the ne- 
glected right part of (f26f would lead to further contribu- 
tions which scale with higher powers of A. 

In the rest of this Section we intend to further simplify 
the above equation and also link the results to those which 
were found in Ref. [TO] . Plugging (l2"2"j) into J2"7|) leads to 



n-l 
3=0 



»>3 



(31) 



where g^r) are the diagonal elements of the matrix 



9 t *(t,h) = ~[v(t),[v(t 1 ),p lJl : 



(32) 



that is, g l a(r) = (i, fx \ g^it, ii) | /x, i). (In the following the 
fixed index /i is suppressed.) We directly obtain 

n— 2 Ti—l n— 2 

Y i C h {r)C k {r 1 ) = H^runr:.^"'-"'-' ■ (33) 

fe=0 i,j=0 fc=0 



where G l (t) is the integral corresponding to (? 4 (t). Since 
Pq(ti) > Pq{t)i we eventually end up with the "best case", 
if we simply set p q (t\) = p q (t) in (f30|) . that is, 



S(t)> 



Jo 4=0 



R(t)R(h). 



(35) 



e. g. 



The "worst case" results by setting p q (t\) = p q (0), 

Pq(tl)/p q (t) < e M 2.7 for i < Tfl. 

It is worth to mention that the above equation does 
not depend on the additional variables. Moreover, since we 
typically consider relaxation times which are much larger 
than the time scale at which correlations decay, G l (t\) and 
R{t\) are approximately constant rates, at least as long as 
correlations do not reappear. We hence obtain 



S(t) > t 



J2G\t) 2 -R(t) 2 



(36) 



This result coincides with the fourth order estimation S(t) 
which were derived for investigations in the context of 
relaxation in closed quantum systems, see Ref. |10j . 



5 Range of Validity of the Second Order 

In this Section we are going to quantify the validity range 
of the second order prediction which is obtained from the 
original projection onto the local densities only. To this 
end we will define a measure \ which is suitable for any sit- 
uation. But in the context of completely decaying and not 
reappearing correlations this measure directly determines 
the range of length scales on which diffusive behavior is 
to be expected, that is, 



'In 



2VX, 



(37) 



where q min (l max ), g max (imin) correspond to the longest, 
respectively shortest exponentially relaxing p q (t). Below 
^min ballistic transport occurs, whereas beyond / ma x any 
non-diffusive type of transport, e. g., insulating behavior 
may emerge. 

To start with, we consider the two contributions which 
occur in ([3H]) . Their ratio 



f(t) 



W 2 S(t) =w s(t) 



WR{t) 



R(t) 



(38) 
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typically is a monotonically increasing function, cf. (|36|) . 
As a consequence there always exists a time ts with f{ts) = 
1, that is, a time where both contributions are equally 
large. But this fact does not restrict the validity of the sec- 
ond order prediction, if ts 3> tr and therefore / (tr) <C 1. 
The validity obviously breaks down only in the case of, say, 
/(tr) w 1 or even larger. 

Since /(£) and tr depend on W, we use the definition 
of the relaxation time 



exp 



W I dti R(ti) 
lo 



exp 



- 1 



(39) 



in order to replace W in (|38[) . Due to this replacement 
I{ t r) becomes a function of the free variable tr, 



f(r R ) = S(r R ) R(r R ) / dhR(h 



(40) 



(tr still depends on W, of course.) Because also /(tr) usu- 
ally turns out to increase monotonically, we define max(rR) 
as the maximum tr for which /(tr) is still smaller than 1. 
Note that this maximum relaxation time already specifies 
the validity range of the second order prediction. 

However, we usually deal with decaying correlations 
and it is hence useful to set max(rR) in relation to Tc- 
We therefore define the measure \ as the dimensionless 
quantity \ = Tc7/max(TR). For example, x = 1 directly 
implies the breakdown of the second order prediction on 
relatively short time scales on the order of tc, whereas 
X = strongly indicates its unrestricted validity. 

For practical purposes an interpretation of x m the 
context of length scales certainly is advantageous. Such 
an interpretation essentially requires the inversion of (|39p . 
In general this can only be done by numerics. But if cor- 
relations decay completely and do not reappear, we have 
tr = 1/(W R) for tr 3> tc- For sufficiently small q we 
may approximate W w q 2 X 2 . Therefore, for fixed A, we 
may write 



f 



1 



W max R 
1 

W~R 



= 2r c , 
_ max(TR) 



X 2 R 



(41) 



where the factors 2 and 1/2 are chosen to slightly fulfill 
tc <C tr <C max(rR), that is, q mul , q ma , x correspond to the 
longest, respectively shortest exponentially relaxing p q {t). 
We finally end up with <7min/9max ~ 2 ^Jx- (Analogously, 
for fixed q, one obtains A m i n /A max w 2 v/x.) 

It remains to clarify what estimation for S(t) should 
be chosen for the calculation of x- This choice basically 
depends on the intention: One may show that x is small 
even in the "worst case" or that x is large in spite of the 
"best case" assumption, see Sec. [H 



5e 
AE 



A A 



H = 



fj, = N-l 



Fig. 1. A chain of N identical, weakly coupled subunits which 
feature a non-degenerate ground state, an energy gap AE and 
an energy band Se with n equidistant states. The dots indicate 
excitation probabilities and are supposed to visualize a state 
from the investigated "single-excitation space" . 



6 Application to Models 

6.1 Modular Quantum Systems with Random 
Interactions 

In the present Section we firstly consider a model which 
is "designed" for the application of our method, since it 
perfectly fulfills almost all properties which have been as- 
sumed in the previous Secs.[3][5] The projection onto local 
densities and second order TCL remarkably allow for a 
complete characterization of all available types of trans- 
port and their dependence on the considered length scale, 
too. 

According to Fig.[TJ the model is a chain of N identical 
subunits which are assumed to feature a non-degenerate 
ground state, a large energy gap AE and a comparatively 
narrow energy band Se with n equidistant states. In the 
following we focus on the "single-excitation subspace", 
that is, only one subunit is excited to its band, while all 
other subunits are in their ground states. Consequently, 
the local Hamiltonian is given by ([2]) with /i-independent 
coefficients h l = iSe/n + AE. 

The next-neighbor interaction ([3]) is also supposed to 
be identical for all adjacent subunits. In particular the 
//-independent coefficients v l ' J form a normalized matrix 
whose elements are chosen at random from a Gaussian 
distribution with zero mean. But we do not intend to apply 
random matrix theory or discuss quantum chaos. For our 
purposes the crucial point is the fact that matrices of this 
kind satisfy the van Hove property, see Sec. [4] 

However, following the ideas of quantum chaos, we ex- 
pect that each sufficiently complex one-particle system 
will take on a form which is similar to our model, once 
it is partitioned into local subunits and those subunits 
are diagonalized. The local interactions will effectively be 
random (Gaussian orthogonal/unitary ensemble) and the 
local spectra will be more or less equidistant (Wigner level 
statistics). The influence of the spectral details is discussed 
below. 

Our system may be viewed as a simplified model for, 
e. g., a chain of coupled atoms, molecules, quantum dots, 
etc. In this case the hopping of the excitation from one 
subunit to another corresponds to transport of energy, es- 
pecially if AE ^S> Se. This system may also be viewed as a 
tight-binding model for particles on a lattice. In this case 
the hopping corresponds to transport of particles. There 
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are n "orbitals" per lattice site but apparently no particle- 
particle interaction in the sense of the Hubbard model. 
Although the are random, these are systems without 
disorder in the sense of, say, Anderson |14J . since the v % 'i 
are independent from /i. For some literature on this class 
of systems, see Refs. [13I15I16I7I17] . Remarkably, a very 
similar model has been used in Ref. [18] in order to in- 
vestigate the flow of wave energy in reverberation room 
acoustics. 

In the context of this model we are mainly interested 
in the probability p^t) for finding an excitation of the 
/xth subunit to its band, while all other subunits are in 
their ground state. This quantity is described by (|3} with 
p l = 1/y/n. Because the model is translational invariant 
and the overall conservation of probability is naturally 
provided, the projection onto these local quantities and 
TCL2 immediately leads to (|20|) . that is, the relaxation 
of p q (t) with the standard rate — WR(t). The underlying 
correlation function C(r) reads 



C(r) 



£ 

i,j=0 



v^l 2 COS 



Se (i - j) t 



(42) 



Of course, C(r) depends on the concrete realization 
of the random coefficients v %, K But, due to the law of 
large numbers, the crucial features are nevertheless the 
same for the overwhelming majority of all realizations, 
at least as long as ^/n ^> 1. And in fact, C(t) typically 
assumes the form in Fig. O It decays like a standard corre- 
lation function on a time scale in the order of tq = 1/Se. 
The area under this initial peak is approximately given 
by R = 2irn/5e. However, because the local bands pos- 
sess an equidistant level spacing, C(r) is a strictly peri- 
odic function with the period T — 2-Kn/Se, unlike stan- 
dard correlation functions. As a consequence its time in- 
tegral R(t) nearly represents a step function, see Fig. [21 
A non-equidistant level spacing will certainly "smooth" 
these steps and change their width, height and distance 
as well. But we expect that the tendency of an increasing 
rate R(t) will nevertheless remain. 




Fig. 2. Sketch of the correlation function C(r) [dashed line] 
and its integral R(t) [continuous line]: C(r) features complete 
revivals at multiples of T such that R(t) has a step-like form. 
This sketch indeed reflects the numerical results for C(r) and 
R(t) but highlights the relevant time scales for clearness. 



Along the lines of Sec. El for tq < t < T diffusive 
behavior with a diffusion constant D = X 2 R is to be 
expected, cf. ([21~j) . And indeed, for density profiles p q {t) 
which decay on such an intermediate time scale we find 
an excellent agreement between (|2ip and the numerical 
solution of the full time-dependent Schrodinger equation 
(which is obtained by incorporating Bloch's theorem and 
exactly diagonalizing the Hamiltonian within decoupled 
subspaces) . In Fig. [3] a "typical" example is shown for a 
single realization of the random numbers v l J . 




0.08 



Fig. 3. Time evolution of a mode p-n(t) which decays on an 
intermediate time scale tc < t <T. Numerics (crosses) shows 
an exponential decay which indicates diffusive behavior and 
is in accord with the theoretical prediction (continuous line). 
Parameters: N = 120, n = 500, 5e = 0.5, A = 0.0005. 



However, until now the above picture is not complete 
for two reasons. The first reason is that p q (t) may decay 
on a time scale that is long compared to T. According to 
([2TJ]) . this will happen, if 2 (1-cosg) X 2 RT > 1 is violated. 
If we approximate 2 (1 — cosq) w q 2 — 4tt 2 / I 2 for rather 
small q (large I), we obtain the condition 



47r 2 n X 
I6e 



> 1. 



(43) 



If this condition is satisfied for the largest possible I, i. e., 
for I = N, the system exhibits diffusive behavior for all 
modes. If, however, the system is large enough to allow 
for some / that violates condition (j43|), diffusive behavior 
breaks down in the long- wavelength limit. This result is 
again backed up by numerics, see Fig. |U 

Towards what transport type does the system deviate 
from diffusive, if condition (|4"3"1) is violated? As already 
mentioned above, we have to consider time scales f > T 
in this regime. We may thus approximate R(t) w 2Rt/T, 
see Fig. [2] Plugging D(t) = 2 X 2 Rt/T into © leads to a 
spatial variance a (t) = 2X 2 Rt 2 /T, clearly indicating a 
transition towards ballistic transport. The validity of our 
approach is again backed up by numerics: In the ballistic 
regime a Gaussian decay of p q (t) is to be expected, see 
Fig. [Hi. 

In a second case transport may be non-diffusive, if the 
p q {t) decay on a time scale that is short compared to tq- 
This will happen, if 2 (1 — cos q) X 2 R tc <C 1 is violated. If 
we use the approximation 2 (1 — cosq) w 4 for the largest 
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Fig. 4. Deviations of the time evolution of Fourier modes p q (t) 
with q — 2n/N, the longest wavelength, from a purely expo- 
nential decay for different model parameters N and n. These 
deviations are based on a measure used in [16] and are in accord 
with the claim that diffusive transport behavior is restricted to 
the regime defined by condition ((43} ■ Other model parameters: 
Se = 0.5, A = 0.0005. 




k/(N-l) 



k/(N-l) 



k diagram for N = oo 
In the limit of small N 



Fig. 6. Upper section of the Ek-vs. 
(continuous lines) and N = 10 (dots) 
the band structure of distinct smooth lines apparently breaks 
down towards a disconnected set of points. Other parameters: 
n = 500, Se = 0.5, A = 0.0005. 




11 0.0013 

t/T t/T 

Fig. 5. (a) Evolution of a mode p^/6o(~t) which decays on a 
time scale t S> T, parameters: N — 120, n = 500, Se = 0.5, 
A = 0.0005. (b) Evolution of a mode Pn(t) which decays on a 
time scale t <C tc, parameters: N = 120, n — 500, Se = 0.5, 
A = 0.004. In both cases numerics (crosses) shows a Gaussian 
decay which indicates ballistic behavior and is in accord with 
the theoretical predictions (continuous lines). 



particle interactions. Nevertheless, in the limit of many 
bands (many orbitals per site) and few sites (few k- values) 
two features may occur: Firstly, the band structure in k- 
space becomes a disconnected set of points rather than the 
usual set of distinct smooth lines, see Fig. [6l It is there- 
fore impossible to extract velocities by taking derivatives 
of dispersion relations. And secondly, the eigenstates of 
the current operator no longer coincide with the Bloch 
eigenstates of the Hamiltonian such that the current be- 
comes a non-conserved quantity, even in the absence of 
impurity scattering. It is straightforward manner to check 
that both features occur in the regime where condition 
(|43|) is fulfilled. This is the regime where standard solid 
state theory breaks down due to the fact that the system 
is too "small" . 



Validity of the TCL Approach and Failure of the NZ Tech- 
nique 



possible q (smallest possible I), the above inequality may 
be written as 



If this inequality is violated, diffusive behavior breaks down 
in the limit of short-wavelength modes. Moreover, if the 
second order still yields reasonable results for not too large 
A, we expect a linearly increasing rate R(t) and thus a 
Gaussian decay, that is, according to the above reasoning, 
ballistic transport. For increasing wavelength, however, 
the corresponding inequality will eventually be satisfied, 
hence allowing for diffusive behavior. Also these conclu- 
sions are in accord with numerics, see Fig. [5)d. 

Relationship to Standard Solid State Theory 

Standard solid state theory always predicts ballistic trans- 
port for a translational invariant model without particle- 



Lowest order TCL suggests the emergence of ballistic trans- 
port, if either condition (|4"3"|) or (|4"4"|) is violated, e. g., if 
the coupling strength A is sufficiently weak or strong, re- 
spectively. Although the agreement of these predictions 
with numerics surely is evident, it is desirable to support 
them by an independent analytical calculation, too. Such 
calculations are indeed possible in the limit of A — > and 
A — > oo, since then first order perturbation theory can 
be invoked: A V is a small perturbation to Ho in the first 
case and vice versa in the second case. In both limits a 
localized initial state ^(0) = eventually leads to a 

variance <x 2 (i) oc t 2 in the limit N — > oo. The correspond- 
ing proof essentially requires the use of Bessel functions 
and their properties, see Appendix lAl for details. 

However, in a sense the agreement of the numerical 
simulations with the TCL2 result is really surprising: The 
fact that the correlation function C(t) features full re- 
vivals at multiples of T points towards strong memory ef- 
fects, cf. Fig. [21 It appears to be a widespread belief that 
long memory times have to be treated by means of the 
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NZ projection operator technique. Whereas the solutions 
of NZ2 and TCL2 are almost identical for t# < T in the 
diffusive regime, for tr T in the deep ballistic regime 
the NZ2 equation contrary predicts a purely oscillating 
behavior of the corresponding density profiles p q (t). But 
such a behavior obviously contradicts the observed Gaus- 
sian decay in Fig. [5] and consequently demonstrates the 
failure of NZ2 in the description of the long-time dynam- 
ics. 




1 2 3 4 5 

tr/T 

Fig. 7. Ratio f(t) of fourth to second order contribution at 
t = tr ("best case" approximation). The fact that /(tr) ~ 1 
is already reached for rather small tr ~ T seems to indicate 
the breakdown of the TCL2 prediction for ballistic behavior 
on large length scales. Parameters: n = 500, Se = 0.5. 



Unfortunately, for this model the unrestricted validity 
of lowest order TCL cannot be prognosticated by the ideas 
of Secs.|4]and[5l According to Fig. at tr ~ T the fourth 
order takes on the same order of magnitude as the second 
order. This finding wrongly indicates the breakdown of 
the TCL2 prediction for tr ^> T. One is apparently con- 
cerned with a situation where higher order contributions 
are not individually small but otherwise compensate each 
other to approximately zero. In such a situation a large 
but finite fourth order term alone is not suitable as an 
"alarm" criterion. Note that for this model the fourth or- 
der term does not increase arbitrarily and converges to a 
finite value, see Fig. [7] This will not be the case in the 
following Sec. 16.21 



6.2 Anderson Model 

Since it had been suggested by P. W. Anderson, the An- 
derson model served as a paradigm for transport in disor- 
dered systems |14ll9l20l21l22l23j . In its probably simplest 
form the Hamiltonian may be written as 

H = e r at a r + ^2 & t *V . ( 45 ) 

r NN 

where the a r , a J are the usual annihilation, respectively 
creation operators; r labels the sites of a d-dimensional 
lattice; and NN indicates a sum over nearest neighbors. 
The e r are independent random numbers, e. g., Gaussian 
distributed numbers with mean (e r ) = and variance 



(e r e r ') = 5 r ^'(J 2 . Thus, the first sum in (|45ll describes 
a random on-site potential and hence disorder. 

It is well known that in the presence of disorder, a 7^ 0, 
the eigenstates of the Hamiltonian (f45|) are no longer given 
by Bloch functions: The eigenstates are not necessarily 
extended over the whole lattice and can become localized 
in configuration space, i. e., the envelope of a wavefunction 
decays exponentially on a finite localization length. 

This phenomenon and its impact on transport have in- 
tensively been studied for the Anderson model [14119120121] 
For the lower dimensional cases, d — 1 and d = 2, it 
is commonly assumed that all eigenstates feature finite 
localization lengths for arbitrary (non-zero) values of a. 
Consequently, in the thermodynamic limit, i. e., with re- 
spect to the infinite length scale, an insulator is to be 
expected, see Ref. |20|21) . Of particular interest is the 3- 
dimensional case, of course. Here, the mobility edge ar- 
ranges spatially localized and extended wavefunctions into 
separated regimes in energy space. When the amount of 
disorder is increased, the mobility edge goes above the 
Fermi level and a metal-to-insulator transition is induced 
at zero temperature [20 21 , still on the infinite length 
scale. When a is further increased, above some critical 
disorder <jc, all eigenstates become localized and an insu- 
lator is to be expected for T > also, (ac — Wc /\/l2 w 6, 
where Wc is the critical "width" of the Gaussian distri- 
bution [2T].) 

However, with respect to finite length scales the follow- 
ing transport types are generally expected: (i.) ballistic on 
a scale below some, say, mean free path; (ii.) possibly dif- 
fusive on a scale above this mean free path but below the 
localization length; and (iii.) insulating on a scale above 
the localization length. 

In the present Section, other than most of the perti- 
nent literature, we do not focus on the mere existence of 
a finite localization length. Instead we rather concentrate 
on the size of the intermediate regime and the dynam- 
ics within. We especially demonstrate that there exists a 
length scale regime in which the dynamics is indeed diffu- 
sive and characterized by an energy-independent diffusion 
constant. In principle, this regime could be very large for 
long localization lengths. But the results in this Section 
indicate that it is not. Investigations in this direction (but 
not for d — 3) are also performed in Refs. [23124) . 

Our approach is still based on the form of the TCL 
technique which has been established in Sees. [3] - In 
this form the TCL method is restricted to the limit of 
infinite temperature. This limitation implies that energy 
dependences are not resolved, i. e., our results are to be 
interpreted as results for an overall behavior of all energy 
subspaces. Thus, the regime in which TCL2 holds is char- 
acterized by the fact that the dynamics within is diffusive 
at all energies with a single diffusion constant. 

In principle the formalism is also applicable in the case 
of finite temperatures, but solely in the limit of weak cou- 
pling strengths. In this limit the energy subspaces of the 
full system are directly known from the spectra of the 
local subunits. The contributions of the initial condition 
p q (0) to energy subspaces can therefore be weighted with 
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a Boltzmann factor and may be treated separately from 
each other. But for strong interactions, as it is the typical 
case for the Anderson model, energy subspaces are simply 
not extractable from the uncoupled system and merely the 
case of infinite temperature is accessible. 



of a correlation function 



H U 
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Fig. 8. A 3-dimensional lattice which consists of N layers with 
n x n sites each. Only next-neighbor hoppings are taken into 
account. Constants for intra-layer hoppings are set to 1 (white 
arrows), inter-layer hoppings are specified by another constant 
A (black arrows). 



C(t) 



1 



Tr{ v^fi+iit) 



(46) 



which depends only negligibly on the concrete choice of the 
layer number \x (during some relevant time scale). Simple 
numerics indicates that this assumption is well fulfilled 
(for the values of a which are discussed here), once the 
layer sizes exceed ca. 30 x 30. Exploiting this assumption 
immediately leads to the TCL2 prediction ([20]) . namely, 
the relaxation of Fourier modes p q (t) with the standard 
rate —WR(t). The underlying correlation function is the 
above C(r), of course. 

Direct numerics shows that C(r) again looks like a 
standard correlation function. We thus retain the former 
notation of tq as the correlation time and R as the area 
under the initial peak of C(r). The numerical results also 
indicate that neither tc nor R depend substantially on n 
(at least for n > 30). Consequently, both tq and R are 
essentially functions of a. 



As shown in Fig.[5J we consider a 3-dimensional (cubic) 
lattice consisting of N layers with n x n sites each. For 
technical reasons we use a Hamiltonian which is almost 
identical to with a single exception: All those terms 
corresponding to hoppings between layers are multiplied 
by some constant A. However, for A = 1 the Hamiltonian 
reduces to the standard Anderson Hamiltonian l|15]). 

A "coarse-grained" description in terms of subunits is 
established now: At first we take all those terms of the 
Hamiltonian which only contain the sites of the /zth layer 
in order to form the local Hamiltonian of the subunit 
[i. Thereafter all those terms which contain the sites of 
neighboring layers fi and fi+ 1 are selected in order to form 
the interaction A w^.^+i between adjacent subunits y, and 
[i + 1. Then the total Hamiltonian may be also written in 
the form of as H = Hq + XV. Note that in this form 
the additional parameter A allows for the independent ad- 
justment of the interaction strength. The eigenbasis of Hq 
may be found from the diagonalization of disconnected 
layers. 

By p^ we denote the particle number operator of the 
/xth subunit, i. e., the sum of a\ a T over all r of the /xth 
layer. Since the overall number of particles is conserved, 
YYli^Pn, H] = 0, and no particle-particle interactions are 
taken into account, we choose to restrict the analysis to 
the one-particle subspace. We may therefore implement 
Pn by (j4|) with p 1 — l/n. The corresponding expecta- 
tion value Pfj,(t) is the probability for locating the par- 
ticle somewhere within the /ith subunit. The considera- 
tion of these "coarse-grained" probabilities corresponds to 
the investigation of transport along the direction which is 
perpendicular to the layers, cf. Fig. [H Instead of simply 
characterizing whether or not there is transport at all, we 
analyze the full dynamics of the p^ (t) . 

According to Sec. [3l the application of our method 
requires that two conditions are fulfilled: (i.) the overall 
conservation of probability which is naturally provided; 
and (ii.) the "average" translational invariance in terms 




Fig. 9. Time evolution of modes p q (t) with q = -n (shortest 
wavelength) and q = 2n/N (longest wavelength) for A = 0.24 
(circles) and A = 0.08 (squares). Remaining model parameters: 
n = 30, a — 1, and N = 10. Symbols are obtained numerically 
from exact diagonalization, whereas all curves correspond to 
the theoretical TCL2 predictions. Additional data from a 4th 
order Suzuki- Trotter integrator is shown in (b) for the case of 
A = 1 and N — 42 (triangles). Note that the curves in (b) and 
(c) are identical, since the factor W is the same, see text for 
details. Further note that the transition from (a) towards (d) 
is to be expected for larger iV with a single choice of A. For 
evidence in that direction, see Fig. 1101 
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Fig. 10. Numerical illustration of the transition in Fig. [9] for 
a fixed A. Model parameters: n = 20, a = 1, N = 26, and 
A = 0.24. Smaller layer sizes (n = 20) are chosen to allow 
for exact diagonalization of a system with a larger number of 
layers (N = 26). Note that theoretical curves are not indicated, 
since for n = 20 the "average" translational invariance is not 
well fulfilled yet. 



According to all above findings and the reasoning in 
Sec. El for t > tc diffusive behavior with a diffusion con- 
stant D = A 2 R is to be expected in TCL2, cf. (plj). Due 
to the independence of R from both n and N, the perti- 
nent diffusion constant for arbitrarily large systems may 
be quantitatively inferred from the diagonalization of a 
finite, e. g., "30 x 30 layer". 

In order to check the above theory, we exemplarily 
present some results here. For n = 30 and, e. g., a = 1 
we numerically find tc ~ 10 and R m 2.9. Thus, addi- 
tionally choosing A = 0.24 and considering the longest 
wavelength in a N = 10 system {q — tt/5), we obtain 
W R 0.064. This corresponds to a ratio tr/tc ~ 1.6 
and hence tr > tc which justifies the replacement of (|20l) 
by (f2"Tj) . And indeed, for the dynamics of p q (i) we get an 
excellent agreement of the theoretical prediction based on 
(|2"Tj) with the numerical solution of the full time-dependent 
Schrodinger equation, see Fig. [9}x Note that this solution 
is obtained by the use of exact diagonalization. Naturally 
interesting is the "isotropic" case of A = 1 . Keeping a = 1 , 
one has to go to the longest wavelength in a A" = 42 sys- 
tem in order to keep the W of the former example un- 
changed. If our theory applies, the decay curve should be 
the same. This indeed turns out to hold, see Fig.[9jD. Note 
that the integration in this case already requires approx- 
imative numerical integrators like, e. g., Suzuki- Trotter 
decompositions |25j . A numerical integration of systems 



with larger N rapidly becomes unfeasible but an analysis 
based on ([21]) may always be performed. 




Fig. 11. Numerical results for the measure \ with respect 
to the disorder a and the inverse layer size 1/n 2 ("best case 
approximation" ) . \ detects the corridor of length scales where 
the overall behavior of all energy regimes is diffusive. For those 
values of \ which are in the order of 1/4 the corridor does 
not exist. But for those values of \ which are closer to 0, the 
corridor opens and diffusive behavior is present in all energy 
regimes with a single diffusion constant. The smaller \, the 
larger is this corridor of diffusive length scales. An absolute 
minimum Xmin ~ 0.02 is found at a ~ 0.5 in the limit of 
n — > oo. Note that only 10% of the whole area is extrapolated 
(the area in front of the thick line) . 

So far, we have characterized the dynamics of the dif- 
fusive regime. We now turn towards an investigation of 
its size. To this end we consider the measure x which has 
been introduced in Sees. H] and O Recall that this measure 

yields l min /l max ~ 2vx, cf. m- 

Imin is the length scale where diffusive dynamics breaks 
down due to the fact that the corresponding p q (t) decays 
on a comparatively short time scale tr ss tq- In com- 
plete analogy to the former model, this transition is still 
correctly described by TCL2: R(t) is not approximately 
constant but linearly increases during the relaxation pe- 
riod. This strongly indicates a transition towards ballistic 
behavior, see Fig. [9K 

Contrary, Z max is the length scale where diffusive dy- 
namics breaks down, because the corresponding p q (t) de- 
cays on a long time scale tr at which the additional con- 
tribution S(t) becomes non- negligible, cf. (|3"0)) . Since the 
interaction perfectly fulfills the van Hove structure, S(t) 
essentially reflects fourth order effects which account for, 
e. g., insulating behavior, cf. Fig. [5JI. [Note that all our 
data available from exact diagonalization is in accord with 
a description based on l[30|) . Note further that S(t), other 
than R(t), scales significantly with n. This eventually gives 
rise to the n-dependence in Fig. ITTI ] 

As shown in Fig. 1111 for each layer size n there is some 
disorder which minimizes x an( f hence "optimizes" the 
diffusive regime. However, for n = 30 (back of Fig. ITTj) 
we find ^Xmin ~ 1/3 at this optimum disorder, indicat- 
ing about one diffusive wavelength. Exactly those respec- 
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tive wavelengths have been chosen for the examples in 
Figs. EbjC but not in Fig. [9jl. For all a and up to n = 100 
(which is about the limit for our numerics) x clearly ap- 
pears to be of the form x(<r, n) = A(a)/n 2 + B{a). Ex- 
trapolating this l/n 2 -behavior yields a suggestion for the 
infinite model n — oo (front of Fig. [TTj) . According to 
this suggestion, we find y/Xmin ~ 1/7, again at optimum 
disorder. This indicates a rather small regime of diffusive 
wavelengths, even for the infinite system. 

We finally recall that these findings apply at infinite 
temperature, i. e., the small diffusive regime is charac- 
terized by the fact that the dynamics within is diffusive 
at all energies with a single diffusion constant. The ab- 
sence of diffusive dynamics in the limit of strong disorder 
apparently agrees with the common expectation that dif- 
fusion at each energy stops completely, once a reaches a 
value larger than ac ~ 6 and all eigenstates become lo- 
calized 20 21j. Down to intermediate amounts of disorder 
the localization phenomenon certainly plays a crucial role 
for the smallness of the diffusive regime. In the limit of 
weak disorder, however, localization appears to be less im- 
portant, especially since only few eigenstates in the outer 
tails of the density of states are localized, while the over- 
whelming majority is extended. In that limit the smallness 
of the regime in which TCL2 holds may result from the 
fact that diffusive dynamics is present at each energy but 
with a energy-dependent diffusion constant, cf. [21]. Then 
the overall behavior of all energy subspaccs is a multi- 
exponential decay and not diffusive, of course. On that 
account it appears plausible that the diffusive regime be- 
comes much larger, if only a certain subspace of energy is 
considered. 



7 Summary and Conclusion 

In this work the TCL projection operator technique has 
been applied to quantum transport in modular systems. 
In particular, the projection onto local densities and the 
lowest order contribution of the TCL expansion have been 
used as a strategy for the analysis of transport and its 
length scale dependence. Furthermore, an estimation for 
the range of validity of the lowest order prediction has 
been derived by means of an extension of the standard pro- 
jection that includes additional relevant degrees of free- 
dom. This estimation is especially suited for interaction 
types with van Hove structure and also provides an inter- 
pretation in the context of length scales. 

As a first application a single-particle model has been 
investigated which has randomly structured interactions 
but nevertheless is translational invariant. For this model 
a full characterization of all types of transport behavior 
has been obtained. Remarkably, purely diffusive behavior 
has been demonstrated to occur on an intermediate length 
scale which is bounded by completely ballistic regimes in 
the limit of short and long length scales. 

The next step in the context of single-particle models 
has been performed by the introduction of disorder. The 
3D Anderson model has been shown to exhibit fully dif- 
fusive dynamics on an intermediate length scale between 



mean free path and localization length. But it has also 
been demonstrated that the diffusive regime is extremely 
small, at least in the considered limit of high tempera- 
tures. 

Naturally, a further important step is given by the ap- 
plication of the theory to interacting many-particle sys- 
tems with or without disorder. In fact, there are already 
preliminary results for a "many-particle" model which is 
translational invariant, has randomly structured interac- 
tions but does not allow for any single-particle restriction. 
This model has been discussed only briefly |16ll5j and ba- 
sically seems to show the same dynamical behavior as its 
single-particle analog, that is, diffusive dynamics on inter- 
mediate length scales breaks down towards ballistic trans- 
port in the limit of large length scales. The main problem 
for this system arises from the fact that a direct compari- 
son of the second order prediction with the numerical so- 
lution of the time-dependent Schrodingcr equation is not 
feasible for sufficiently large systems. Consequently, an ap- 
propriate estimate for the validity range of second-order 
TCL of the type derived in this paper is indispensable. 

A Analytically Solvable Limit Cases 

The results of Sec. 16.11 strongly suggest the occurrence of 
ballistic transport in the two limit cases (|43|) and (|44|) 
of weak and strong coupling strength A. Since these re- 
sults depend on the applicability of the TCL2 method, it 
would be desirable to support them by independent ana- 
lytical calculations. Such calculations are possible in the 
two limit cases A — > and A — > oo, since then first order 
perturbation theory can be invoked. Moreover, we will re- 
strict ourselves to the case where the Hamiltonian given 
by Eqs. (Jl])-([3|) is translational invariant and the next- 
neighbor interaction matrix (|3|) is symmetric, i. e., we will 
assume 

^ = hl = h\ (47) 

<m+i=<i=0^'. (48) 

In these calculations various sums over fj, = 0, 1, . . . , N — 1 
occur, which will be most conveniently approximated by 
integrals. This approximation is exact in the limit N — > 
oo. The variance ([6]) will be considered with local "exci- 
tation densities" §4§ using p 1 ^ = 1. 

Let us assume that the symmetric next-neighbor interac- 
tion matrix v has been diagonalized: 

n-i 

5>^ = Sfe S* (49) 

We will consider the time evolution of an initially localized 
excitation, i. e., a solution 

¥>(*) =£>!*(*) lAM> (50) 

of the Schrodinger equation with the initial value ip l (0) = 
S/j.,0 fiifi- 
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References 



In the limit A — > and iV -> oo we obtain 



2 ~S h0 \J M (2\v°>°t)\ 2 



(51) 



where Jui denotes the |/x|-th Bessel function, and as a 
consequence 

CT 2 (*)=5> 2 |^)r~2A 2 ,W, (52) 
indicating ballistic transport. 
A. 0.2 The case A — > oo 

Upon rescaling the Hamiltonian |T|) in the form 

H=\h + V (53) 
A 

we can again apply first order perturbation theory. This 
yields eigenvalues of the form 



where 



1 2?r/i 
Ek,n = j E k + 2 s fc cos 



E k = J2\Sl\ 2 h i 



(54) 
(55) 



Consequently 



YM S l exp(^|^ ]./|„|i 2 s,/) 



(56) 



In order to evaluate the variance ([6]) we will utilize the 
following formula, stated without proof, 



^ m 2 J n+m (s) J m (t) 



171= — OO 

+2 



= -j [J n +2(s - t) + J„_ 2 (s - i) + 2 J„(s - t)] 
+ 1 [-4+1 (s - i) - J n -i (« - *)] 



(57) 



in the special case of n = 0. We then obtain for the leading 
term 

a 2 (t)~2t 2 Y / \S°k\ 2 4 + 0(t), (58) 

k 

again indicating ballistic transport. Numerical tests show 
that the parabolic approximation (l58|) is very close for all 
times. 
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